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ABSTRACT 

Invisible plasma content in blazar jets such as protons and/or thermal electron- 
positron (e^) pairs is explored through combined arguments of dynamical and ra- 
diative processes. By comparing physical quantities required by the internal shock 
model with those obtained through the observed broadband spectra for Mrk 421, we 
obtain that the ratio of the Lorentz factors of a pair of cold shells resides in about 
2 ^ 20, which implies that the shocks are at most mildly relativistic. Using the ob- 
tained Lorentz factors, the total mass density p in the shocked shells is investigated. 
The upper limit of p is obtained from the condition that thermal bremsstrahlung 
emission should not exceed the observed 7-ray luminosity, whilst the lower limit is 
constrained from the condition that the energy density of non-thermal electrons is 
smaller than that of the total plasma. Then we find p is 10^-10'^ times heavier than 
that of non-thermal electrons for pure pairs, while 10^-10^ times heavier for pure 
electron-proton (e/p) content, implying the existence of a large amount of invisible 
plasma. The origin of the continuous blazar sequence is shortly discussed and we spec- 
ulate that the total mass density and/or the blending ratio of pairs and e/p plasma 
could be new key quantities for the origin of the sequence. 
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1 INTRODUCTION 

The discovery of strong inverse Compton components in X 
and 7-ray emission from jets in active galactic nuclei (here- 
after AGN) for a wide range of spatial scales (e.g., Coll- 
mar 2001 for review) enables us to probe quantitatively 
the energetics of relativistic jets. The kinetic power of non- 
thermal electrons has been estimated by various authors 
both for inner core jets (i.e., blazars) (e.g.. Kino, Takahara 
and Kusunose 2002, hereafter KTK; Kusunose, Takahara 
and Kato 2003) and large scale jets (e.g., Tavecchio et al. 
2000; Leahy and Gizani 2001, 2002; Kataoka et al. 2003). 
However, the material content of relativistic jets is not easily 
constrained by observations since the emission is dominated 
by that from non-thermal electrons and probably positrons 
and it is difficult to directly constrain thermal matter con- 
tent. Hence, the plasma composition in AGN jets, whether 
normal proton-electron (e/p) plasma or electron-positron 
pairs (e*) is a dominant composition, is still a matter of 
open issue (e.g., Reynolds et al. 1996; Celotti, Kuncic, Rees 



and Wardle 1998; Wardle et al. 1999; Hirotani et al. 1999; 
Sikora and Madejski 2000; Ruszkowski and Begelman 2002; 
Kino and Takahara 2004, hereafter KT04). This problem 
prevents us from estimating the total mass and energy flux 
ejected from a central engine. 

To constrain invisible matter content such as thermal 
electron-positron pairs and/or protons co-existing with non- 
thermal electrons, dynamical considerations are indispens- 
able. In KT04, we proposed a new procedure to constrain 
the invisible thermal plasma component in classical FR II 
radio sources. We used the fact that the mass and energy 
densities of the sum of thermal and non-thermal particles 
are larger than those of non-thermal electrons which are de- 
termined by observations. Here we apply the same technique 
to the inner core jets of AGNs (i.e., blazars) based on the in- 
ternal shock model. The internal shock model is believed to 
be most plausible to explain the production of high energy 
photons and time variabilities in blazars. It has been widely 
applied also to the prompt emission of gamma-ray bursts 
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(hereafter GRBs) (e.g., Rees 1978; Rees and Meszaros 1994; 
Kobayashi, Piran and Sari 1997;, Daigne and Mochkovitch 
1998; Ghisellini 1999; Spada, Ghisellini, Lazzati and Celotti 
2001). It is worth to note that recently GhiseUini et aL (2005) 
proposed a structured jet model cosisting of a fast spine 
surrounded by a slowly moving layer for explaining VLBI 
scale radio blobs. At the present, however, it is not evident 
where is the acceleration site of electrons in the structured 
jet model. This is one of the prime issues which should be 
answered. Internal shocks are potentially the building blocks 
of the spine part of the structured jet. Whereas we recog- 
nize the importance of the detailed structure of jets, as a first 
step we focus on the physical condition of the flow based on 
the simple internal shock model. 

The methodology of constraining the invisible plasma 
content in the emission region is as follows. As mentioned 
above, a lower limit to the total mass density (sum of non- 
thermal electrons and invisible plasma) is restricted by the 
definition that the mass density of total plasma should be 
smaller than that of the non-thermal electrons. The mass 
density of non-thermal electrons can be estimated by multi- 
frequency observations. For this purpose, in ij3]we review the 
shock dynamics of two colliding shells. Note that we do not 
use the simple two point-mass approximation (e.g., Piran 
1999; Lazzati et al. 1999; Zhang and Meszaros 2004 for re- 
view) but employ the exact shock dynamics throughout this 
work. This makes outcomes more accurate. In SjJJwe briefly 
review the previous results on the amount of non-thermal 
electrons based on KTK. In ^ we constrain on the amount 
of total mass density. As for the upper limit, we use the con- 
straint that bremsstrahlung emission from thermal electron 
(and positron) component should not exceed the observed 7- 
ray emission. We postulate synchrotron self-Compton (SSC) 
emission dominance in the 7-ray band which is supported 
by the observed correlations between TeV7-ray and X-ray 
in TeV blazars (e.g., Takahashi et al. 1996, 2000; Catanese 
et al. 1997; Maraschi et al. 1999). We can thus bracket the 
amount of total mass density in the emission region from 
below and above. In this way we apply this method to the 
archetypal TeV blazar Mrk 421. In SjB] we further estimate 
the shock dissipation rate of the colliding cold shells. The 
dissipation rate is a widely discussed quantity in literatures 
concerning gamma-ray bursts (e.g., Lazzati, Ghisellini and 
Celotti 1999; Piran 1999). The shock dissipation is believed 
to be the ultimate source of heating and accelerating parti- 
cles. Summary and discussion are in ^ 



2 KEY FEATURES OF THIS WORK 

The key features of this work are briefly summarised here in 
advance. The existence of copious amount of invisible plasma 
is predicted from a qualitative consideration. 

2.1 Existence of invisible plasma content 

As mentioned in the introduction, we constrain on the 
amount of invisible plasma content by introducing the dy- 
namical considerations. The point is that we divide total 
mass and energy densities into two components, i. e., those 
of non-thermal electrons and those of the other invisible 
components. The comparison of these obtained quantities 



enables us to constrain on the amount of invisible plasma 
and this is a new attempt compared with the previous works. 

Bearing this in mind, next we show a quantitative con- 
sideration which derives the existence of invisible plasma in 
colliding shells of blazar jets. Let us discuss a collision be- 
tween a pair of equal mass-density shells for instance. In 
the comoving frame of one shell, the particles of the other 
shell are coming in with a relative bulk Lorentz factor Tij 
(see in §3 for details) of a few at most (shown in Table 1). 
When only pair plasma are present and all of them are 
accelerated, then the average Lorentz factor of non-thermal 
electrons (je) is expected to be {7e) ~ Pij. This is too small 
to account for observed blazar spectra {'je) ~ 300 (this is the 
case of Mrk 421) obtained by KTK. Therefore only a frac- 
tion of the pair should be accelerated, the ratio of the 
rest mass density of non-thermal electrons to that of total 
plasma is about ^ij/i'Je)- Similarly we can discuss the case 
for shells with pure e/p plasma makeup. If all of the dis- 
sipated energy goes into the electron acceleration, then we 
have {7e) ~ [mp/m^)Vij. This is too large to account for the 
spectra and it requires a limited fraction of electrons being 
accelerated. Thus, invisible plasma is qualitatively expected 
when the internal shock is responsible for the production of 
non-thermal electrons. In this paper, we will quantitatively 
explore the amount of invisible plasma in jets. 

2.2 Why we use shock dynamics? 

In the previous studies, the colliding shells have been ap- 
proximately modeled as the simple two-point-mass collision 
(e.g., Piran 1999; Lazzati et al. 1999; Zhang and Meszaros 
2004). The reason why we use the shock dynamics instead 
of the two-point-mass model is as follows. When one try to 
derive the mass density from the mass, one eventually needs 
to know lengths and velocities and they can be consistently 
obtained by the shock model. Hence the shock analysis is 
the best way for investigating the invisible plasma content 
in jets. 



3 SHOCKS IN COLLIDING SHELLS 

Here we review the relativistic shock jump conditions. We 
use one-dimensional shock dynamics of a pair of colliding 
shells to apply the standard internal shock model to blazars. 
Suppose the situation in which a rapid shell overtakes a pre- 
viously ejected slow shell. There are four characteristic re- 
gions designated by (1) unshocked slow shell, (2) shocked 
slow shell, (3) shocked rapid shell, and (4) unshocked rapid 
shell. These regions are separated by the forward shock (FS), 
the contact discontinuity (CD), and the reverse shock (RS). 
In this paper, we use the terminology of regions i (i=l, 2, 
3, and 4) and position of discontinuity i (i=FS, CD, and 
RS) where FS, CD, and RS stand for the forward shock 
front, contact discontinuity, and reverse shock front, respec- 
tively. The fluid velocity and Lorentz factor in the region i 
measured in the interstellar medium (hereafter ISM) frame 
are expressed as Vi(= fiic) and Ti, respectively. The rela- 
tive velocity and Lorentz factor of the fluid i measured in 
the frame j are denoted by Vij{= —Vji — (iijC = —Pjic) 
and Fji), respectively. Rest mass density, pressure, 

and internal energy density are expressed as pi, Pi, and ei, 
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respectively. As for the equation of state (EOS), we take 
Pi = (71: — ^ PiC^)^ where 7^ is the adiabatic index. We 
sometimes use the subscripts s and r instead of 1 and 4, such 
as Ti = and r4 = T^. 

In the hmit of strong shock, with the assumption of 
cold upstream (Pi = 0), the jump conditions for the forward 
shock are written as follows (Blandford & McKee 1976): 



p2 _ (ri2 + i)[72(ri2-i) + i]' 

J- FSl — 



62 = ri2P2 



72(2-72)(ri2-l) + 2 ' 

P2 _ 72ri2 + 1 

pi 72 — 1 ' 



(1) 



where ri2 = rir2(l — (iiP2), and Fpsi is the Lorentz factor 
of forward shock measured in the rest frame of the unshocked 
slow shell. In the relativistic limit, the adiabatic index is 
72 = 4/3. Using the same assumptions as in the forward 
shock, the jump conditions for the reverse shock are given 
by: 

(r34 + i)[73(r34-i) + i]' 



p2 

J- RS4 



73(2 - 73)(r34- l)+2 ' 

pi l-A^-iA + 1 

63 = f 34/53 , = — ^ 

pA 73-1 



(2) 



where = r3r4(l — (3^/34), and rRS4 is the Lorentz fac- 
tor of the reverse shock measured in the rest frame of the 
unshocked rapid shell. The equality of pressure and velocity 
across the contact discontinuity gives 



P2 =P3 



(3) 



After the shocks break out the shells, r2 = is not 
satisfied because a rarefaction wave changes the density and 
velocity profiles (e.g.. Kino, Mizuta and Yamada 2004, here- 
after KMY). We do not treat the rarefaction waves for sim- 
plicity, concentrating on the major duration before shock 
breakout. It may be useful to rewrite the pressure balance 
along the CD as 

£4 ^ (72ri2 -n)(ri2 - 1) 
pi (73r34 + i)(r34-i) ■ 



(4) 



In general, the number of physical quantities in each 
region is 3, pi, Pi (or e^), and Vi. Forward and reverse shock 
speeds (i.e., vfs and «fs) are two other quantities. In all, 
there are 3 x 4 + 2 = 14 physical quantities. Note that Pi and 
a are connected with EOS. The total number of the jump 
conditions is 3-1-3-1-2 = 8. Hence, given 3-1-3 = 6 upstream 
quantities for each shock, we can obtain the remaining 8 
downstream quantities by using 8 jump conditions. It is to 
be noted that the absolute value of the rest mass density is 
irrelevant to the shock dynamics since the shock dynamics 
is linear with respect to the mass desnity. Then, actually 
we need to specify 5 quantities if we give the density ratio 

P4/pl. 

For a specific case for TeV blazars, we here impose the 
following two conditions; (i) the unshocked shells are cold, 
i.e., Pi = P4 = 0, (ii) the Lorentz factor of the shocked 
regions F3(= F2) is identified as that of the emission region 
obtained by the observed broadband spectra. Further, we 
examine the following three cases for the ratio pr/ps', (a) the 
energy of bulk motion of the rapid shell {E = Fmc'^) equals 
to that of the slow one in the ISM frame (we refer to it as 
"equal energy (or E) case"), (b) the mass of the rapid shell 
(m = pFA) equals to that of the slow one (hereafter we call 



it "equal mass (or m) case"), and (c) the rest mass density 
of the rapid shell equals to that of the slow one (hereafter 
we call it "equal rest mass density (or p) case"). Here, A 
denotes the thickness of the shell measured in ISM frame. 
These choices are based on the conjecture that the ejecta 
from the "central engine" is likely to have a correlation with 
each other (e.g., NP02; KMY). Hereafter, we assume that 
the widths of two shells are the same in the ISM frame, that 
is Ar/As = 1 (e.g., NP02, Spada et al. 2001). Note that in 
the case of Ar = As and Fr > Fs, Ps is always larger than 
pr for equal E and equal m cases. 

Thus, we give 4 quantities, Pi, P4, F2 = F3 and one 
relation between the rapid and slow shells, p^j p\ depending 
on cases (a) through (c) described above. As a remaining 
quantity, the Lorentz factor of the rapid shell F4 is treated 
as a free parameter. Although we do not specify the absolute 
value of p, we treat the abosolute value in actual applica- 
tions. It is compared with that of non-thermal electrons in 
the shocked regions as described in !j5l The absolute value of 
the rest mass density comes into play when two-body pro- 
cesses such as bremsstrahlung emission is used to obtain the 
upper limit of p. We will properly discuss these points. 

In the following sections, we focus on the values of (i) 
the value of Fi and F4, (ii) 63 and/or ps, as a tool to examine 
the physical quantities of invisible matter content. 



4 AMOUNT OF NON-THERMAL ELECTRONS 

4.1 Number and energy densities 

Based on the detection of inverse Compton emission in 7-ray 
band, the number and energy densities of the non-thermal 
(hereafter "NT") electrons r^'^ and e^"'^ in shocked regions 
can be determined by the comparison of the observed broad- 
band spectrum and the theoretical one. Although the mini- 
mum Lorentz factor of relativistic electrons is not definitely 
determined and affects mainly the number density n^^, 
we regard that low energy electrons below 7e,min consti- 
tute thermal electrons. Considering the observed fiat num- 
ber spectrum of electrons, fixing 7e,min = 10 does not cause 
any major problem with n^^. 

Here, we briefiy quote the resultant n^^ and e^^ ob- 
tained in KTK. Hereafter, we omit the subscript express- 
ing the regions i(= 2, 3) for simplicity. For clearness of the 
following argument, we define that n^^ and e^^ also in- 
clude NT positrons when they exist. The quantity n^^ is 
written as ri^'^ = ne(7e)d7e, while e^''" is given by 

<^'^ = {7e)n^^meC^, where ne(7e) and {7e) are the energy 
spectrum and the average Lorentz factor of NT electrons, re- 
spectively. By a detailed comparison of the SSC model with 
observed broadband spectrum of Mrk 421, we obtained 

as 

-0.6 



NT 



11 X 



cm 



(5) 



Here, we adopt the index of injected electrons for Mrk 421 
as s = 1.6 (e.g., Mastichiadis & Kirk 1997; Kirk & Duffy 
1999) and the case of 7e,min = 10 was examined in KTK. 
The best choice of the size of the emission region is P = 
2.8 X lO^^cm with an order of magnitude uncertainty. Thus, 
the corresponding uncertainty of n^"^ amounts to two orders 
of magnitude; for smaller P, larger n^""" should be adopted. 
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But, as far as the the shock dynamics is concerned, only the 
density ratio plays a role, therefore we adopt the above value 
as the canonical one. 

As for the average energy of NT electrons, we obtained 



NT / NT 



: 3.1 X IQ^m^c. 



(6) 



Since for a = 1.6, electrons near the cooling break energy 
~ 7e,br carry most part of the kinetic energy and e^"^ has 
a weak dependence on 7e,min provided that 7e,min is smaller 
than 7e,br ~ 10'*. Note that the case of 7e,min ~ 10"* is ruled 
out for Mrk 421 since the case does not fit the EGRET data 
(KTK). Therefore, Eq. ([6]) is justified in any case for Mrk 
421. 



4.2 Forward and reverse shocks 

To clarify whether the observed non-thermal emission comes 
mainly from FS or from RS region, the typical frequency 
of non-thermal synchrotron radiation and internal energy 
density in each region are examined here. 

According to the standard diffusive shock acceleration, 
the acceleration time scale is estimated as (e.g., Drury 1983) 
iacc = (27r7emec5)/(e_B) where ^ = A/r^ is a parameter 
related to the amplitude of magnetic fluctuations, A and 
rg are the mean free path for the scattering of electrons 
and Larmor radius, respectively. Here the shock speed is 
taken to be c. On the other hand, the synchrotron cooling 
time is given by tayn = {Q-Km^c?) / {(jT'yecB'^) . The maximum 
Lorentz factor of the non-thermal electrons is evaluated as 



7max oc B ^^'^ by using the condition of 



teyn at 7 

n 



with the assumption that ^ in FS and RS regions takes the 
same value. Hence, the characteristic synchrotron photon 
energy is given by /iJ/syn,o,max oc TiBy^,^ = const., Hence, 
the value i^syn.o.max in FS region and RS regions is the same. 

The total internal energy of NT electrons in FS and 
RS regions may be discussed as follows. If 3> 1 and 

^ 1 are satisfied, we have €2 + P2 — 63 + -P3. In the 
actual case of blazars, Fij is close to order of unity and we 
have P2 ~ Pa- Thus, the energy densities of regions 2 and 
3 are similar and the internal energy is controlled by the 
comoving shell widths. Since r43 ^ is always satisfied, 
co-moving length of RS region is larger than that of FS 
region in the case of Ar = Ag (e.g., Kobayashi and Sari 
2001; NP02; KMY). Thus, the radiation from RS dominates 
over that from FS region. Based on this consideration, we 
focus on RS dominated case in this paper. Hereafter, we 
omit the subscript 3 for simplicity. 



5 CONSTRAINTS ON THE AMOUNT OF 
INVISIBLE PLASMA 

5.1 Lorentz factors of cold shells 

It is hard to estimate the bulk Lorentz factors of cold shells 
simply because they are invisible. However, by using the 
value of Lorentz factor of shocked shell which corresponds 
to the beaming factor of the emission region, we can con- 
strain on the Lorentz factors of the cold shells. Here, we 
consider the range from 3 to 100 for Fr and Fg. Following 
Begelman, Rees & Sikora (1994), we consider the upper limit 
of the Lorentz factors of the emission region as Fr.max = 100, 



Table 1. Lorentz factors obtained by internal shock analysis for 
Mrk 421 



case 


F, 


F, 


Fia 


r43 


equal p (largest Fr/rs) 


3 


48.0 


2.125 


2.125 


equal p (smallest Fr/Fs) 


8.485 


16.97 


1.060 


1.060 


equal m (largest Fr/Ps) 


5.12 


100 


1.35 


4.22 


equal m (smallest Fr/Fs) 


8.983 


17.959 


1.042 


1.082 


equal E (largest Fj/Fs) 


8.57 


100 


1.057 


4.226 


equal E (smallest Fr/Fg) 


9.48 


18.96 


1.027 


1.106 



Notes: r„ 



100, r„ 



■■ 3, Fr/Fs > 2, and r2 = F3 = 12, are 



employed in this analysis. 



while as for the lower limit we employ Fg^min = 3 based on 
Wardle & Aaron (1997). Here, we exclude cases of very weak 
collisions with Fr/Fs < 2 as in NP02. As for the adiabatic 
index in Eq. Q, we approximate 7^ = 4/3 for Tij > 2, oth- 
erwise 7i = 5/3 for simplicity (e.g.. Kirk and Duffy 1999). 

For the TeV blazar Mrk 421, we have already obtained 
F2 = Fs = 12 by the observed multi-frequency spectrum 
(KTK). Hence, Eq. Q is solvable for Fs given Pr/ps and F^. 
Qualitatively, a faster Fr requires a slower Fs to attain the 
same value of F3. Thus, minimun value of Fr corresponds to 
Fr/Fs — 2, while the maximum value of Fr corresponds to 
Fr = 100 or Fs = 3. 

In Table [1] we show the minimum and maximum values 
of Fr and Fs and the corresponding relative Lorentz factors 
F12 and F34 which control the shock heating of the down- 
streams (see Eqs. ((T} and ((2])). From this, we see that the 
range of Fij lies between 1.03 and 4.2. In other words, a 
mildly relativistic shock is realized in the case of Mrk 421. 
We also note that our adopted value of 7e,min = 10 is a rea- 
sonable choice with the assumption that 7e,min should be a 
few times larger than Tij. The corresponding value of Fr/Fs 
is found as 



2 < Fr/Fs < 16.0 
2 < Fr/Fs < 19.5 
2 < Fr/Fs < 11.7 

respectively. 



(equal p), 
(equal m), and 
(equal E), 



(7) 



5.2 Total mass density 

5.2.1 Lower limit of p 



In 15.11 we show that shock is at most mildly relativistic 
though each shell moves at a relativistic speed. As a conse- 
quence, dissipation efficiency is relatively small and (7e) ^ 
F34 is satisfied. Therefore e^^/e = (7e)p^^/F34p < 1 gives 
a tighter constraint than p^^/p < 1. By rewriting the con- 
dition of e^^/e < 1, the lower limit of the total mass density 
is given by 



P 

_NT 



(Te) 
F34 



3.1 X 10^ 



(8) 



Here we omit the subscript of region number i = 3 for 
the various densities for thumbnail writing. From this we 
directly see that in order to accelerate electrons up to 
{■Je) ~ 3.1 X 10^ in the framework of standard internal shock 
model, where only a small available shock dissipation energy 
F34 ~ a few is realized, the invisible mass density at least 
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about 100 times the rest mass density of NT electrons is 
definitely required. In other words, we need a loading of 
baryons and/or a thermal pair plasma. It is worth to note 
the effects of an uncertainty with {7e). We estimated the 
uncertainty range as 2.3 x 10^ < {je) < 4.3 x 10^ (KTK). 
The uncertainty simply leads to a shift of lower limit curve 
by the same factor. Since it causes only a small change on 
the resultant value, we focus on the best-fit case in this work 
for simplicity. 



electron-electron bremsstrahlung is not considered in these 
equations. It is clear that the upper limit of p in this case is 
mp/me times larger than n^nie. 

Lastly, let us check the timescale of pair annihi- 
lation tann. It is evaluated as tann — &e/{neO-TC) ~ 6 X 

cm ^) ^ sec. Hence we see that the annihila- 
tion time scale is much longer than the dynamical time scale 
Uyn = VSR/c « 2 X 10'^(i?/10^'* cm) sec. Therefore pair 
annihilation is not effective in this situation. 



5. 2. 2 Upper limit of n J 

Here, we constrain the upper limit of the number density of 
thermal electrons nj. As mentioned in the Introduction, it is 
widely accepted that observed GeV and TeV 7-rays are SSC 
dominated. In the MeV range, bremsstrahlung radiation by 
the thermal electrons with temperature Oe = kTe/me(? ~ 

~ a few MeV is expected if adequate amount of thermal 
electrons exist in the emission region. At the moment, we do 
not have any observational evidence for the bremsstrahlung 
in MeV band. At the same time, it is fair to note that ob- 
servation in MeV range itself is a challenging area (e.g., 
Takahashi et al. 2003). Here we estimate the upper limit 
of the number density of thermal electrons by assuming the 
observed bolometric luminosity of bremsstrahlung Lbrom.o 
should be lower than that of SSC Lssco which is estimated 
as Lssc.o = 7 X 10** erg s"' (KTK). 

For e^e~ plasma content, we employ Eqs. (21) and (22) 
of Svensson (1982) which express the emissivity of relativis- 
tic e^e~ bremsstrahlung ebrcm(0e, wj) where nj is the num- 
ber density of thermal electrons. Note that these expressions 
do not include the bremsstrahlung between electron-electron 
and positron-positron and the limit will be severer by a 
factor of ~ 2 if we include them. Then, the condition of 
issc.o > ibrom,o is rewritten as 

nj < 9.7 X 10^ [6^^(1-1- 1.76^'')]"'''" cm'^ (Oe < 1) 

< 5.7 X 10^ [ee(ln(l.iee) -^5/4]"^/^ cm"^ (Gc ^ 1). 

(9) 

The bolometric luminosity of the optically-thin 
bremsstrahlung is estimated by Lbrem.o = (47rJ?^/3)r3ebrom 
with the emission size 7? = 2.8 x 10^^ cm and the Lorentz 
factor Fs = 12 as obtained by the broadband spectral fitting 
of Mrk 421 (KTK). The electron temperature is evaluated 
by (73 — l)0e = r34 — 1. The upper limit turns out to be 
about a thousand times larger than the mass density of 
non-thermal electrons. It is consistent with and relatively 
close to the required lower limit of the mass density by Eq. 
([§}. This upper limit depends on the adopted value of R, 
and it is proportional to R~^^^ . Considering that is 
roughly proportional to R~^ , the ratio of this upper limit 
to only has a weak dependence on R. 

Similarly, in the case of electron-proton (hereafter e/p) 
plasma content, we can rewrite the condition of Lssc.o > 

-^brcm,o aS 

nj < 9.5 X 10^ [ey^(l -I- 1.786'-^*)] '^''^ cm"^(ee < 1) 

< 9.6 X 10^ [ee(ln(l.iee + 0.42) + 3/2] ^^^^ cm"^ 

(Se ^ 1) (10) 

with Eqs. (17) and (18) of Svensson (1982). Note that 



5.2.3 Allowed range of p 

We thus obtained the upper and lower limits on p/p^"^ and 
the results are shown in the plane of mass density of in- 
visible plasma and Fr/Fg in the cases of "equal p", "equal 
m" , and "equal E" in Figs. [1] [J] and O respectively. They 
are obtained by solving Eq. ^ and inserting F34 into Eqs. 
([HI, ((9]), and (|10p . The qualitative features are the same for 
these three cases, although different in quantitative detail. 
Summing up in advance, the most important result is that a 
large amount of mass density of invisible plasma is required 
in the emission region. As the value of Fr/Fs increases, the 
value of F34 becomes larger and the lower limit on the in- 
visible mass density {p/p^'^) reduces. Below we discuss two 
extreme cases of different plasma content. One is the case of 
the jet with pure pair plasma content, whilst the other 
is the jet made of pure e/p plasma. 

For pure pair jet, the resultant total mass density 
normalized by is 

2 X 10^ < p/p^'^ < 2 X 10^ (equal p) 
7 X 10^ < p/p^'^ < 2 X 10^ (equal m) 

6 X 10^ < p/p^'^ < 2 X 10^ (equal E). (11) 

For the jets consisting of pure plasma, the predicted 
p/Pe'^ is constrained in a narrow range around 100-1000 
as shown in Figs. [T] [5] and O The number density fractions 
of the shock accelerated pairs are directly obtained as 
Pe'^/p = n^'^/(n<F + n^'^) ~ 10"^ - 10"^ This seems a 
reasonable result since the number of accelerated particles 
is expected to be a small fraction of the thermal pool. 

In the case of pure e/p content, the allowed range of 
p/p^'^ are found to be 

2 X 10^ < p/p^'^ < 3 X 10^ (equal p), 

7 X 10^ < p/p^'^ < 3 X lO** (equal m), and 

6 X 10^ < p/p^'^ < 3 X lO** (equal E), (12) 

respectively. The maximum values of p/p^"^ are about 
nip /me times larger than those in the case of pure pair 
content. 

5.3 Allowed range of e/e^^ 

As shown above, the lower and upper limit of p/p^'^ have 
been obtained in 15.2.11 and 15.2.21 respectively. By using the 
obtained p/p^"^ shown in 15.2.31 we can estimate e/e^'^ = 
F34p/{7e)p^^. For the case of pure content, since the 
allowed range of p/p^"^ is narrow, the corresponding e/e^^ 
is also well constrained as 

1 < e/e^^ < 7 (equal p), 

1 < c/e^^ < 7 (equal m), and 



6 M. Kino and F. Takahara 



1 < e/e^"^ < 7 (equal E), 



(13) 



where we employ Eq. and Table 1. Thus we find that 
e/e^'^ < 2 X 10^ X 1.1/310 ^ 7. In other words, for pair 
content, the total kinetic power of the shocked (emission) 
region is less than Ljjin ~ 7Lj^;^ ^ where ^ is the kinetic 
power of NT electrons estimated as L^^^^^ = 4 x 10''* erg s~' 
(KTK). In the case of e/e^'^ « 1, the non-linear dynamical 
structure of the shock (e.g., Drury and Voelk 1981; Berezhko 
and Ellison 1999) is required for analysing the phenomena at 
the vicinity of the shock front. Note that the case discussed 
here is consistent with our choice of Tr/Ts > 2. 

On the contrary, for pure e/p content, the energetics 
relevant to thermal electrons and NT and thermal protons 
is all quite uncertain. Based on Eq. p2|l and Table 1, we can 
derive 



1 < e/er < 1 X 10" 
1 < e/e^^ < 1 X 10'' 
1 < e/e^^ < 1 X 10' 



(equal p) 
(equal m) 
(equal E). 



(14) 



For the case of maximum values of p/p^"^ in Eq. (|12|) . the 
total kinetic power for pure e/p content reaches Lkin ~ 
10''l/JJ'in g which is extremely large and unlikely. 



6 ON THE SHOCK DISSIPATION RATE 

In order to examine the allocation of the bulk kinetic energy 
of cold shells -Ebik into the thermal energy, we estimate the 
shock dissipation rate of bulk kinetic energy of colliding cold 
shells. Here we denote the thermal energy of shocked shells 
as. E — Eq where Eq is the rest mass energy and E is the 
total kinetic energy which satisfies E oc V2e2 + Faea. Then 
the shock dissipation rate ediss defined as the ratio of the 
thermal energy of mass elements after the collision to that 
of bulk kinetic energy of mass elements before the collision 
is given by 



^diss — 



V2{E-Eq) 



£blk — 
F2KF12 



' l)5m2 + (F34 - l)<5m3] 



(Fi - l)<5m2 + (F4 - l)5m3 



(15) 



where 5m2 and 5mz are the surface mass of the shocked re- 
gions 2 and 3, respectively. Here, 5m,2 and Sm^, are expressed 
as 5m2 — F2p2(wFS — vcD)St and Sma — F3p3(iicD — VRs)St 
where 5t is the corresponding duration time in the ISM 
frame. Two difi'erences between the present work and the 
two-point-mass collision model (e. g., Piran 1999) are that 

(i) we estimate ediss with the shock junction conditions, and 

(ii) we subtract the irreducible rest mass term from the de- 
nominator. Our definition is superior to the previous one 
when the value of relative Lorentz factor (i.e., F34 and/or 
F12) are close to order unity and/or a small Lorentz factors 
for cold shells. From Eq. (|15p we obtain 



0.07 < Ediss < 0-44 (equal p) 
0.09 < Ediss < 0.63 (equal m) 
0.06 < Ediss < 0.35 (equal E). 



(16) 



The larger (smaller) Fr/Fa becomes, the larger (smaller) ediss 
realizes. As previously mentioned (Kobayashi and Sari 2001; 
KMY), the case for equal mass of colliding shells realizes 



largest value of ediss and an asymmetry of each mass reduces 
the value of ediss- 

The fraction of 1 — ediss of the bulk kinetic energy of the 
cold shells survives and transferred to a larger scale. This is 
responsible for the large scale structure such as radio lobes 
and cocoons. Therefore the comparison with the large scale 
kinetic power such as an extended radio emissions of blazars 
(e.g., Antonucci & Ulvestad 1985) will be an important fu- 
ture work although this is beyond the scope of this work. 



7 SUMMARY AND DISCUSSION 

Invisible plasma content in blazar jets such as protons 
and/or thermal pairs is investigated. In this work, we di- 
vide total mass and energy densities into two components, i. 
e., those of non-thermal electrons and those of the other in- 
visible components. It enable us to constrain on the amount 
of invisible plasma in the jet. This is a significant forward 
step compared with previous studies. 

The methodology of constraining the invisible plasma 
content in the emission region is as follows. The lower limit 
to the mass and energy densities of total plasma is limited 
by the definition that the mass and energy densities of total 
plasma should be larger than those of the non-thermal elec- 
trons. The total mass and energy densities are constrained 
by the internal shock dynamics. On the other hand, the up- 
per limit of mass and energy densities for non-thermal elec- 
trons are constrained by the condition that bremsstrahlung 
emission from thermal electron (and positron) component 
should not exceed the observed SSC 7-ray emission. We can 
thus bracket the amount of total mass and energy densities 
in the emission region from below and above. 

We apply this method to the archetypal TeV blazar Mrk 
421 and obtain the following results. 

(1) Mildly relativistic shock is realized. 

By imposing the condition of the bulk Lorentz fac- 
tor of the emission region as F3 = 12 estimated by the 
multi- frequency spectrum of Mrk 421 (KTK), we explore 
the allowed range of Fr/Fs within the framework of the 
standard internal shock model. Adopting the conditions of 
Fs.min > 3, Fr.max < 100, and Fr/Fs > 2 based on the 
literatures (Wardle and Aarons 1997; Begelman, Rees and 
Sikora 1994; NP02), we find that the values of Fr/Fs for Mrk 
421 are limited in the ranges of 2 < Fr/Fs < 16 (equal p), 
2 < Fr/Fs < 19.5 (equal m), and 2 < Fr/Fs < 11.7 (equal 
E), respectively. As mentioned in Kirk and Duffy (1999), a 
very hard injection index of s ~ 1.6 observed in Mrk 421 
well agrees with this mildly relativistic shock regime (See 
Fig. 3 in their paper). Hence we conclude that mildly rel- 
ativistic shocks take place in Mrk 421 from the analysis of 
the observed spectrum and the internal shock dynamics. 

(2) The mass density of invisible plasma is much heavier 
than that of non-thermal electrons. 

Using the condition that the mass and energy densi- 
ties of non-thermal electrons should be lower than those 
of the total ones, we derive the lower limit of total mass 
density at the shocked region. Since the relative Lorentz 
factor between the shocked and unshocked regions is ex- 
pected to be a few (in Table 1), copious amount of mass 
density of invisible plasma is inevitably required. The upper 
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limit of is constrained by the condition that the luminos- 
ity of bremsstrahlung emission should be smaller than the 
observed 7-ray luminosity which is well explained by the 
synchrotron-self-Compton emission. Combining them, the 
allowed ranges of p/p^'^ for pure pair content arc found 
as 2 X 10^ < p/pe'^ < 2 X 10^ (equal p), 7 x 10^ < p/p^"^ < 
2 X 10* (equal m), and 6 x 10^ < p/p^'^ < 2 x 10* (equal E), 
respectively. For pure e/p plasma content, the upper limit 
of p/pe'^ turns out to be 3 x 10®. 

Although the specific index s — 1.6 for Mrk 421 is dis- 
cussed here, we emphasize that the value s < 2 is common 
character for TeV blazars as they indeed display the smaller 
s than 2 (Kirk and Duffy 1999 for review). For instance, 
the choice of s — 2.2 leads to the synchrotron omission with 
vFi, oc Such a soft spectrum significantly conflicts with 
the observed synchrotron emission in blazars (e. g., Fossati 
et al. 1998). 

(3) Electron acceleration efficiency in the shocked region 

is evaluated. 

Once p/p^'^ is bounded as shown in Figs. 1, 2, and 3, 
we can obtain the electron acceleration efficiency as e/e^^ = 
r34p/(7e)p^^ for given r34. Since the allowed value of p/pnt 
is in the narrow range for the case of pure content, wc 
obtain the electron acceleration efficiency as 1 < e/e^^ < 7. 
Correspondingly the total kinetic power of the emission re- 
gion Lkin resides in the range 1 < Lkin/ikil.e < 7. It is clear 
that a loading of proton component significantly enlarges 
Lkin. For e/p plasma content, the value of Lkin is evaluated 
as 1 < Lkin/Lkin,e < W'^ ■ the case of this mgiximal Lkin, 
too large Lkin could lead to a problem for the energy, al- 
though we do not have a tight constraint on the amount of 
proton loading. 

(4) The shock dissipation rate of bulk kinetic energy of 
colliding cold shells is examined. 

In §5, we further estimate the dissipation rate of bulk 
kinetic energy of colliding cold shells into the internal one via 
the shocks. It is qualitatively clear that the larger (smaller) 
Fr/Fs becomes, the larger (smaller) ediss realizes. The re- 
sultant shock dissipation rate for the colliding cold shells 
resides in the range eaiaa ~ 6-60%. 

Lastly, let us discuss the origin of blazar sequence (Fos- 
sati et al. 1998) which is tightly connected to the nature of 
the central engine. In Fossati et al. (1998) they computed 
average spectral energy distributions from radio to gamma- 
rays for complete sample of blazars. The resultant spectra 
show a continuity in that (i) the synchrotron peak occurs in 
different frequency for different samples/luminosity classes, 
with most luminous blazars peaking at lower frequencies; 
(ii) the peak frequency of the gamma-ray component corre- 
lates with the peak frequency of the lower energy one; (iii) 
the luminosity ratio between the high and low energy com- 
ponents increases with bolometric luminosity. They claimed 
that the continuous sequence of properties may be controlled 
by a single parameter, related to the bolometric luminosity. 
Below we enlighten another new ingredient for the origin of 
sequence. Flat spectrum radio quasars (FSRQs) have the or- 
der of magnitude larger kinetic power and energy dcnsitity of 
the external radiation field than BL Lacs (e.g., Sikora et al. 
1997). Hence, the leptonic components in FSRQs ejecta un- 
dergo stronger radiation drag effect in larger external radia^ 
tion fields (Sikora and Wilson 1981; Phinney 1982; see also 



Iwamoto and Takahara 2002). However, the bulk Lorentz 
factors in FSRQs are comparable to or even slightly larger 
than the ones in TeV blazars in spite of being subject to 
much stronger radiation drag (e.g., Kubo et al. 1998; Spada 
et al. 2001; Kusunose et al. 2003). In order to realize larger 
kinetic powers and larger bulk Lorentz factors against the 
strong radiation drag, we may take a new conjecture that a 
larger baxyon loading may occur for FSRQs. Summing up, 
not only the strength of the external radiation field but also 
the total amount and/or blending ratio of pair and e/p 
could be new key quantities to explore the origin of the con- 
tinuous blazar sequence. 
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10^ : lower limit 
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Figure 1. The allowed region of the amount of mass density of 
total plasma normalized by that of non-tlicrmal electrons p/p^"^ . 
The gray region is allowed for pure pair plasma for "equal p" case. 
The line denoted by e/p shows the upper limit for e/p plasma. 
Horizontal axis shows the ratio of the Lorentz factor of a rapid 
shell to a slow one which lies in 2 < Fr/Fs < 16.0. 
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Figure 2. The allowed region of the amount of mass density of 
total plasma normalized by that of non-thermal electrons p/p^"^ . 
The gray region is allowed for pur pair plasma for "equal m" case. 
The line denoted by e/p shows the upper limit for e/p plasma. 
Horizontal axis shows the ratio of the Lorentz factor of a rapid 
shell to a slow one which lies in 2 < Fr/Fs < 19.5. 
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Figure 3. The allowed region of the amount of mass density of 
total plasma normalized by that of non-thermal electrons p/p^'^. 
the gray region is allowed for pure pair plasma for "equal E" case. 
The line denoted by e/p shows the upper limit for e/p plasma. 
Horizontal axis shows the ratio of the Lorentz factor of a rapid 
shell to a slow one which lies in 2 < Fr/Fs < 11.7. 



